Crystallization of an exciton superfluid 
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Indirect excitons - pairs of electrons and holes spatially separated in semiconductor bilayers or quantum 
wells - are known to undergo Bose-Einstein condensation and to form a quantum fluid. Here we show that 
this superfluid may crystallize upon compression. However, further compression results in quantum melting 
back to a superfluid. This unusual behavior is explained by the effective interaction potential between indirect 
excitons which strongly deviates from a dipole potential at small distances due to many-particle and quantum 
effects. Based on first principle path integral Monte Carlo simulations, we compute the complete phase diagram 
of this system and predict the relevant parameters necessary to experimentally observe exciton crystallization in 
semiconductor quantum wells. 

PACS numbers: 71.35.Lk, 03.75.Hh, O5.30Jp 



I. INTRODUCTION 

Quantum coherence of bosonic particles is one of the most 
striking macroscopic manifestations of the laws of quan- 
tum mechanics governing the microworld. The discovery of 
Bose-Einstein condensation in atomic vapors^ was followed 
by the observation of condensation of bosonic quasiparti- 
cles in condensed matter - excitons. Here we mention early 
claims (though highly controversial) for three-dimensional 
(3D) semiconductors, 2 electron bilayers in a quantizing mag- 
netic field,™ exciton-polaritons in microcavities 5 6 and so- 
called indirect excitons formed from spatially separated elec- 
trons and holesPltLU Not only the bosonic gas phase was ob- 
served but also the formation of a quantum Bose liquid - an 
exciton superfluid with its peculiar loss of friction - could re- 
cently be verified. 4 6 Thus it is tempting to ask whether there 
exists also a solid phase of bosons. 

The key properties of a crystal are particle localization and 
long-range spatial ordering. To achieve spontaneous crystal- 
lization requires to find a Bose system with sufficiently strong 
and long range pair interaction (here we do not consider par- 
ticle localization induced by an external field in an optical lat- 
tice or cavity^^j). However, the vast majority of previous ex- 
perimental investigations have been performed in the regime 
of weak nonideality, where the interaction energy is small 
compared to the quantum kinetic energy. Therefore, promis- 
ing candidates for a bosonic solid are atoms or molecules with 
dipole interactiorp3 or excitons. Here, indirect excitons offer a 
number of attractive features: a strong dipole-type interaction, 
the suppression of biexciton or trion formation, the compara- 
tively long radiative life time (on the order of microseconds) 
and the external controllability of the density and dipole mo- 
ment via an electric field perpendicular to the quantum well 
planepEUm 

In this paper we present clear evidence for the existence 
of a crystal of indirect excitons in semiconductor quantum 
wells. We compute its full phase diagram and reveal the pa- 
rameters for its experimental verification. Our predictions are 
based on first principle path integral Monte Carlo (PIMC) sim- 
ulations. But in contrast to previous quantum Monte Carlo 
studies which predicted crystallization in model systems such 



as el ectro n-hole bilayers, I^ED or two-dimensional dipole sys- 
temsP^EH here we use realistic parameters typical for indirect 
excitons. In particular, we fully take into account the finite 
quantum well width, the composite character of the excitons 
and the different masses of electrons and holes. This turns out 
to be of crucial importance for the exciton-exciton interaction 
which strongly departs from a dipole potential at small dis- 
tances. As a direct consequence we observe that the exciton 
crystal exists only in a finite density interval and undergoes 
quantum melting both at high and low density. Furthermore - 
when the exciton superfluid crystallizes to form a solid, quan- 
tum coherence is lost abruptly, i.e. there is no supersolid ex- 
citon phase. 

This paper is organized as follows. In Sec.|Il]we introduce 
the system of indirect excitons and present its reduced quasi- 



2D description. In Sec. Ill the effective exciton-exciton in- 
teraction potential is derived and its accuracy is verified. In 



Sec. IV we present our simulation results and the phase dia- 
gram of indirect excitons. Finally, we draw our conclusions in 
Sec.M 



II. MODEL 

We consider a semiconductor quantum well (QW) of width 
L containing N e = Nh electrons and holes in the conduction 
and valence band, respectively, which are created by an opti- 
cal pulse. 19 Application of an electrostatic field of strength E 
perpendicular to the QW plane created e.g. by a tip electrode 
allows to spatially separate electrons and holes to different 
edges of the QW. By varying E this separation can be changed 
between and L giving rise to a variable dipole moment d. At 
the same time, the field also provides lateral confinement and 
a variable particle density, via the quantum confined Stark ef- 
fect, for details of the setup, see K. Sperlich et al. 15 Finally, 
the system is kept in thermal equilibrium at a finite temper- 
ature T which does not exceed a few percent of the binding 
energy of an electron-hole pair, thus all electrons and holes 
will be bound in N = iV e indirect excitonsP^ 

The thermodynamic properties of this system are fully de- 
scribed by the density operator of N e electrons and holes, 
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Pn n = A{e~^ H /Z}, where Z is the partition function, 
P = l/(k B T) and A denotes full anti-symmetrization among 
all electronic and hole variables. The full Hamiltonian H 
contains kinetic energy, the interaction with the external elec- 
tric field and all Coulomb pair interactions between the 2N e 
charged particles 



H = H\\ + H z + W 



with the single particle contributions 
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Here denotes inter-particle distances in the QW plane, 
yQW j s tne q W confinement, U is the electrostatic poten- 
tial energy due to the electric field and e is the background 
dielectric constant, m" and m 1 - are the the effective electron 
(hole) masses which take into account the anisotropy of the 
in-plane (out of -plane) parabolic dispersions in the QW. 

Under the present conditions of strongly bound indirect 
excitons with parallel dipole moments resulting in a strong 
exciton-exciton repulsion the very complicated evaluation of 
the density operator p^ N can be substantially simplified. 
As was shown in Ref ~U the system can be mapped onto N 
excitons which can be treated as composite spin polarized 
bosons 22 where deviations from the Bose statistics (arising 
from the original Fermi statistics of electrons and holes) have 
been found negligible. 9 Thus, the density operator is reduced 
to a fully symmetric one of N excitons, p s N . Furthermore, 
all pair interactions can be properly averaged along the QW 
width giving rise to an effective (d-dependent) exciton-exciton 
interaction Vxx- As a result the system 2D Hamiltonian enter- 
ing p s N becomes 
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where to-x = + m|j is the in-plane effective mass, r, the 
in-plane center of mass (com) coordinate of the i th exciton 
and Tij = \ri — Tj\ denotes the com distance between two 
excitons. 



III. EFFECTIVE INTER-EXCITON INTERACTION 

To verify the approximation |4]) and the validity of the po- 
tential Vxx we consider the two-exciton (biexciton) problem. 



We define the exciton interaction energy as the energy dif- 
ference of a biexciton and two single excitons, -Exx(fhii) = 
E2x( r hh)~2Ex, which depends parametrically on the distance 
between the holes in a biexciton problem, = \R\ — R 2 \. 
The distance remains a well defined quantity also at small 
inter-exciton separations, when a strong overlap of the exciton 
wavefunctions and particle exchange takes place. In this case 
the com distance is not physical. The substitution of in 
Eq. Q by r^h can be justified as follows. 

Similar to the hydrogen problem, the single exciton wave 
function can be factorized into the com and the relative part 



*(r,R) = * c (R°)%(\r-R\) 
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where the vectors r, R and R a denote the electron, hole and 
com coordinates, respectively. 

The relative part \& r can be found by solving a single parti- 
cle problem with the reduced mass /i = m^m^j (m e + toj,) 
in the potential, Vd = — e 2 / y/\r — R\ + d 2 , where the z- 
direction is taken into account explicitly by the exciton dipole 
moment d. For the spatially indirect exciton we approximate 
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i.e. the leading term of the expansion describes a harmonic 
oscillator and the relative part near the exciton origin decays 
as a Gaussian. Now, using the definition of Rq and the sub- 
stitution, (r — R) = 7 m (R° — R) with j m — mx/m e , the 
relative part can be expressed solely in terms of the hole coor- 
dinate (keeping the com coordinate R as a fixed parameter) 



*(r, R) = * c (r, R) %(R, R°) 
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where the relative part ^ contains a factor 7^ in the expo- 
nent, tyf ( r )\ r =\R Q -R\ e~ 7 ™ r /' . For a typical electron- 
hole mass ratio in semiconductors, -f m ~ 2 ... 4, we con- 
clude, that the hole is well localized around the com. This 
allows to make a second step. 

We treat the excitons in the Born-Oppenheimer (BO) ap- 
proximation and apply the adiabatic transformation for the 
spatial part of the full wavefunction {Spin degree of free- 
dom are omitted in the present analysis, as this requires a sig- 
nificantly more elaborated simulations. The model used for 
the exciton interaction potential, is assumed to have a signif- 
icantly larger effect on the results, when the spin fluctuations 
in the ferromagnetic phase 2 ^! 

**x = 7^)2 E (^) 5P ' +m %(P,rx,P e r 2 ,R 1 ,R 2 ) 



x^ b (P h Rx,P h R 2 ) , (10) 



which can be symmetric or antisymmetric depending on the 
symmetry of the spin part. The action of the electron and 
hole permutation operators, -P e (h), explore all exchange pos- 
sibilities (excluding the electron-hole exchange). Within this 
ansatz one can self-consistently solve the Schrodinger equa- 
tions for electrons 

H e ^\ ri ,r 2 ,Ri,R2) 

= E^(R 12 )^ n \r 1 ,r 2 ,R 1 ,R 2 ) , (11) 

and holes 



[H h + E^(R 12 )\ ^(Ru^Rl, R° 2 ) 



E { ^^\R^R2,RlRl), (12) 



where 



[t: + V eh (r z - fix) + V eh ( ri - R 2 , 
+ V ee {r 1 - r 2 ) , 

^ = ^ + Hh(J?12) , 

2=1,2 

h 2 v 2 



(13) 
(14) 



e(h) 



2m, 



e(h) 



with n,m 6 { A, S } being defined by the symmetry of the 

(n) 

electron (hole) wavefunction, and Eg ' being an additional 
mean-field electron potential influenced by the holes in the 
biexciton. 

If the holes are treated as infinitely heavyp^the numerical 



solution of Eq. ( 12 1 is not necessary and the biexciton energy 



can be decomposed, Eix = -E-i"^ (?"hh) + V\±, with Vhh = 
e 2 / l^hhl- The electron contribution E^ is the solution for a 
singlet (triplet) state 
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with the holes located at R\. 2 = i^ry,- This equation has 
been solved numerically for an experimentally feasible e-h 
separation d — 13.3agp3A first observation is that, the en- 
ergy E e (m<a — > oo) is not sensitive to r^, once rht ^ d, see 
Fig.[T^. This is understood from the behavior of the electron 
density (see Fig. [TJ)): in all cases the electron cloud extends 
well beyond r^h which is a result of the shallow interaction po- 
tential, V e h(r), of an electron with the two holes for r^h < d, 
and the strong e-e repulsion that keeps the electrons at an av- 
erage distance f ~ 20a% apart, practically independent on the 
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Figure 1: (Color online) Exciton interaction potential Vkx for a 
dipole moment d = 13.3cib. (a) The interaction potential Vxx 
(point-dashed line) is compared to the exciton interaction energy 
Exx in several approximations: average interaction of two exci- 
tons (Exx) evaluated by PIMC simulations using the Hamiltonian, 
Eq. Ml; the BO-model with the infinite hole mass, symmetric (Exx) 
and antisymmetric (Exx) electronic states; the improved BO-model, 
-Exx ~ E e + Vhh, with a realistic mass ratio mh/m e — 2.46 (ZnSe- 
based QW). Also shown are the electronic contribution E e and the 
dipole potential d 2 /r 3 . Two vertical lines indicate the boundaries 
of the exciton crystal, (b) Radial electron density n c (r) for several 
hole separations r^, relative to the mid-point of two holes located at 
(15) .Ri,2 = ±|rhh- (c) Electron pair distribution function g te (r ee ) for 



the values ru, in (b). 



hole-hole separation. This behavior is evident from the pair 
distribution function g{r ee ), see Fig.[T];. Consequently, for a 
large exciton dipole moment, we observe no noticeable differ- 
ence in the energy of the symmetric and antisymmetric states, 
merging into a single curve E e (mh — > oo), see Fig. la. With 
these results we can now analyze -Exx(fhh), cf. red dashed 
line in Fig.[lJ. At large distances, r > d, Exx practically co- 
incides with the classical dipole potential, V D = d 2 /r 3 , so we 
expect the system to behave like 2D polarized dipoles, at low 
densities. At smaller distances, r < d, however, Exx essen- 
tially follows a Coulomb potential which arises mainly from 
the hole-hole repulsion. Finally, for r <C d, the interaction 
energy shows an unphysical Coulomb singularity originating 
from the assumption of an infinite hole mass. In real systems, 
Exx is expected to be softer, approaching a finite value at zero 
distance, due to quantum diffraction and exchange effect, sim- 
ilar to behavior of the Kelbg potential in 3D electron-ion plas- 
masP^EHTherefore, we proceed with the generalization of the 
model for a finite hole mass. 

In the situation with a large dipole moment, as considered in 
Fig.[lJ, the interaction energy is positive at all distances and, 
hence, no bound states (biexcitons) are formed. This origi- 
nates from the positive eigenvalues of the Schrodinger equa- 
tion for the holes ( |12) , Therefore, evaluation of the interaction 
energy should not be limited only to the ground state solution 
of Eq. ( 12 1, but should include contribution of all states, in- 
cluding the continuump2l This can be done directly via the 
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two-particle partition function Z 2 , 
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the density matrix, and the thermodynamic energy estimator 
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E(r hh ) = -— lnZ 2 (/3, r hh ) 
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Here, Z 2 parametrically depends on the distance r\± between 
the particles. Applied to the case of two holes in the biexciton 
(E = E 2 x), the density matrix is the solution of the two-body 



Bloch equation with the Hamiltonian, Hh+E t 
see Eq 
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density matrix and the relative part 
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Here U ef{ is the effective pair actionP^ES introduced in a way 
that at large distances and (or) high temperatures it reduces 
to /3{e 2 / |r hh | + £'e ( " ) (r hh )). Substituted in Eq. (18l-(19|l we 
obtain 
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E 2 x{r hh ;P) = k B T+ {k B T+^U ea (r hh ,r hh ;l3)j , (22) 

where the first term accounts for the com kinetic energy (in 
2D). For spherically symmetric potentials the effective ac- 
tion and its temperature derivative can be evaluated with the 
matrix-squaring technique! 30 ! 31 ! The resulting interaction en- 
ergy, E X x(rhh;P) = ^2x(rhh; P) - E x (0), evaluated at the 
temperature 1//3 = 10~ 3 Ha* is shown in Fig. la by the red 
solid line. Quantum effects arising from the finite hole mass 
(e.g. for the ZnSe-based QWs, mb/m e ~ 2.46) strongly af- 
fect the interaction energy Exx for r < 3a B , which conse- 
quently approaches a finite value at zero distance. 

For final comparison, we compute the exciton interaction 
energy by PIMC simulations using the Hamiltonian Q. We 
used two bosonic excitons of mass m x in periodic boundary 
conditions. The result, (Exx)> as a function of the average 
inter-exciton distance, (r) = J dr rg(r)-(J dr g(r)) _1 , eval- 
uated via the exciton pair distribution function g(r), is shown 
in Fig. la by the solid squares. This quantity agrees well with 
the finite-mass BO solution, Exx, for rj,h > 5a' B , and con- 
firms applicability of both models in the density range where 
we predict formation of the excitonic crystal. The deviations 
being noticeable at smaller distances are outside the density 
range used in the present analysis. 
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Figure 2: (Color online) Constant density freezing. Top row: 2D 
PDF g{rij) [relative to a fixed particle in the center] for naj 2 = 
0.0035 and temperatures fc B T/Ha* of 1.74- 10" 3 (a), 1.38- 10~ 3 (b) 
and 1.08 ■ 10~ 3 (c). Bottom: Radial distribution function \g(r) — 1 
(left) and bond angular order distribution function ga(r) (right) at 
na% 2 = 0.0022 for N = 500. Lines are guide to the eye to visualize 
an algebraic decay in this log-log plot. 



citons can be efficiently computed with full account of all in- 
teractions, quantum and spin effects, without further approxi- 
mations. Below we use atomic units, i.e. lengths will be given 
in units of the electron Bohr radius, ag = h 2 e/(e 2 m'^), and 
energies in units of the electron Hartree, Ha* = e 2 / (ea%). Of 
central importance for the crystallization is the coupling (non- 
ideality) parameter, i.e. the ratio of interaction energy to ki- 
netic energy. For a quantum system with Coulomb (dipole) in- 
teraction it is given by the Brueckner parameter r s (the dipole 
coupling parameter D), 
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where a is the mean inter-particle distance and n the exciton 
density. Note the opposite scaling of r s and D with density. 

We perform 2D grandcanonical PIMC simulations^! with 
periodic boundary conditions and extract the results for the 
canonical ensemble with TV = 60 . . . 500 excitons. To map 
out the phase diagram we scan a broad parameter range span- 
ning three orders of magnitude of density and temperature. 
We first obtain the phase diagram for a fixed value of the 
dipole moment, corresponding to d = 13.3 a B , and after that 
analyze in Sec. [VB] how the crystal phase boundary changes 
when d is varied. 



IV. SIMULATION RESULTS 



A. Spatial ordering of excitons 



Using PIMC simulations with p s N and the Hamiltonian (|4| 
the thermodynamic properties of the TV strongly correlated ex- 



To detect crystallization we compute the exciton pair distri- 
bution function [PDF], g(r). This function is homogeneous in 
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Figure 3: Temperature dependence of the defect fraction at density 
nag 2 = 2.2 ■ 10~ 3 and particle numbers N = 501 - 505 (verti- 
cally aligned dots) characterizing the ordered and disordered phases. 
The defect fraction (1 — Pa) shows a sharp jump at T ~ 10~ 3 Ha*. 
This is in disagreement with the KTNHY theory, which predicts a 
continuous unbinding of dislocations in the hexatic phase indicated 
by a continuous variation of the critical exponent rie(T) < 1/4 and 
the bond angular correlation function ge(r) ~ r ~ r 'e(T) _ j n contrast, 
we observe an abrupt transition from the LR angular order (T < T c ) 
to a quasi-long-range angular order with rje{T) ~ 2 (T > T c ), see 
Fig.0 



an ideal gas, whereas in the fluid and crystal phase it exhibits 
increasing modulations which signal localization and spatial 
ordering. Typical examples of g(r) are displayed in the top 
rows of Figs. [2] and [4] and show clear evidence of exciton lo- 
calization. The existence of the translational long range order 
(LRO) is detected from the asymptotic behavior of the angle- 
averaged function g(r) for large r = \r\. In 2D a possible 
freezing scenario is given by the Kosterlitz-Thouless-Nelson- 
Halperin- Young (KTNHY) theory (see the overview^!), pre- 
dicting an exponential (algebraic) decay of the peak heights 
of g(r) above (below) the melting temperature. Indeed, our 
simulations find some support for this scenario, see bottom 
left part of Fig. |2] 

The existence of angular hexagonal LRO follows from the 
asymptotic behavior of the bond angular correlation function, 
ffe(0 = WS(r)V 8 (0)>. with Ve(rfc) = n^YZi* 



i66 fc 



where m is the number of nearest neighbors of a particle lo- 
cated at j-fc, and Qki is their angular distance. We observe 
a change from an exponential asymptotic of g§ to a con- 
stant which is the expected behavior for a liquid-solid tran- 
sition, see bottom right part of Fig. [2] There are some in- 
dications for the existence of an hexatic phase - coexistence 
of angular quasi-LRO (algebraic decay) and missing transla- 
tional LRO in a narrow temperature interval, see curves for 
k B T = 1.05 • 1(T 3 Ha* and k B T = 1.25 ■ 1CT 3 Ha*. 

In addition we performed a Voronoi analysis, which pro- 
vides access to local distortions of the hexagonal symmetry 
of the lattice. The average fraction of particles (the probabil- 
ity) with a number of nearest neighbors deviating from 6 is 
referred to as the defect fraction, i.e. (1 — P 6 ). The results 
of Fig. [3] explore the nature of the melting transition at con- 
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Figure 4: (Color online) Isothermal freezing and melting of indirect 
excitons. a)-f): 2D PDF g(r) for fc B Ti = 0.001 Ha* at densities 
naf of 0.84- 10" 3 (a), 1.3- 10" 3 (b), 1.7- 10~ 3 (c), 3.2- 10~ 3 (d), 
3.6 ■ 10~ 3 (e), and 4.0 ■ 10~ 3 (f). Bottom panel: Superfluid fraction 
7s, Eq. {15\ , vs. density for two temperatures. Symbols are PIMC 
results, lines are a guide to the eye. The increase of 7 S at high density 
extends over a small finite range of solid-liquid coexistence which is 
due to the finite particle number in the simulations. 



stant density. We observe a sharp increase of the number of 
defects at the melting point which is in disagreement with the 
KTNHY scenario. A possible alternative to the KTNHY is 
a first order solid-liquid phase transition, with an exponential 
decay of g§ (r). However, the latter was not observed in our 
simulations, possibly, due to a limited system size (N ~ 500). 
The constructed Voronoi map for different particle configura- 
tions, shows the accumulation of the defects at the boundaries 
between few crystallites. A similar picture, but for a signifi- 
cantly larger classical system (N ~ 10 6 ) has been recently re- 
ported and the transition was proved to be of the first order.El 
If that system was equilibrated sufficiently long, the interme- 
diate hexatic phase completely vanished. With our data for the 
limited particle numbers we can not give a confident answer 
whether we observe a discontinuous transition in the present 
system. 



B. Exciton quantum coherence. Superfluidity 

After analyzing emergence of spatial ordering let us turn 
to the quantum coherence properties of nonideal indirect ex- 
citons. In a 2D Bose system cooling leads to sudden emer- 
gence of coherence in the liquid phase - the normal fluid - 
superfluid transition. The phase boundary is gov erned by the 
Berezinskii-Kosterlitz-Thouless (KT) scenario 35 and is given 
by the condition x — 4/7s f° r the exciton quantum degener- 
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Figure 5: (a) Temperature dependence of the winding number 
( W 2 ) (T) for the exciton numbers TV = 56 and 170. Density na*, 2 = 
5 ■ 10~ 3 . The Berezinskii-Kosterlitz-Thouless transition tempera- 
ture T KT (N) is determined by the condition, 1 ^ (W 2 )(T K t) = 4,/ir, 
shown by the horizontal dashed line, (b) System size dependence of 
Tkt(L) for three densities: naf = 5 ■ 10" 3 , 10~ 2 and 2 ■ 10~ 2 . 
Values of Tkt are rescaled to fit into a single plot. 



acy parameter x = " A 2 
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where n s = 7 s n is the exciton superfluid density. Therefore, 
a key quantity is the superfluid fraction 7 S , where < 7 S < 1. 
In PIMC simulations, it is directly computed from the statis- 
tics of the winding number 
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Typical simulation results for y s are shown in the bottom part 
of Fig. g| 

Figure [5] illustrates the computation of the winding num- 
ber versus temperature (left) and the finite size scaling for the 
critical temperature Tkt of the BKT transition (right). One 
observes a systematic shift of Tkt (AT) to lower values with 
an increase of the system size N. The extrapolation to the 
thermodynamic limit, T KT (L — > oo), with L = yWV/n, 
is obtained by fitting the simulation data by the equation 
Tct(T) = Tkt(oo) + 6/ln 2 (X) . It is a direct consequence 
of the Kosterlitz-Thouless renormalization group analysis 35 
which is considered to be exact in the asymptotic regime of 
large L. This scaling allows us to make predictions for the 
phase transition line in a macroscopic system. 



C. Phase diagram of indirect excitons 

We now summarize our findings in the complete phase di- 
agram of indirect excitons in the density-temperature plane 
which is presented in Fig. [6] The degeneracy line \ = 1 sepa- 
rates the regions of classical (above the line) and quantum be- 
havior (below). While classical excitons exist only in a fluid 
(or gas) phase the quantum region is composed of three differ- 
ent phases: a normal fluid, a superfluid and a crystal phased 
Correspondingly, there exist two triple points, at the upper left 
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Figure 6: (Color online) Phase diagram of 2D indirect excitons with 
d — 13.3 ag. Circles and squares mark our PIMC results, data 
for triangles are fronP. Vertical dashed lines (D = 17 ± 1 and 
r s = 9.4 ± 0.3) indicate the two density induced quantum freezing 
(melting) transitions. Filled symbols mark the two triple points. The 
normal fiuid-superfluid phase boundary is marked by the red line and 
is below the ideal estimate Tkt according to Eq. {14\ , cf. thick solid 
line labeled \ — 4. The line T& p marks the freezing transition of 
a classical 2D dipole system. The e-h plasma phase is beyond the 
present analysis. 



(right) edge of the crystal phase. At high temperature the ex- 
citons are in the fluid phase. Cooling leads either into the su- 
perfluid or crystal phase. There is no cooling transition from 
the superfluid to the crystal. 

At low densities cooling always leads into the superfluid 
phase; the transition is accompanied by a sudden increase of 
7 S from zero to a finite value. The phase boundary is substan- 
tially below the upper limit Tkt(« s = n), Eq. (24 1, and is in 



full agreement with our previous analysis for 2D dipoles^ 1 in- 
dicating that the exciton interaction is close to a dipole poten- 
tial. The picture suddenly changes when the density exceeds 
na^ 2 ps 0.00078: the superfluid transition vanishes and, in- 
stead, a strong modulation of the PDF is observed signaling 
crystallization, cf. top row of Fig. [2] The critical density cor- 
responds to a dipole coupling parameter D c = 17 ± 1 which 
agrees with studies of pure 2D dipole systems^ 3 !!. Note that 
the freezing temperature changes non-monotonically exhibit- 
ing a maximum value T max around na^ 2 ps 0.002. 

The superfluid-solid transition is verified by simulating 
compression along several isotherms. At low temperature 
and low density, the superfluid fraction 7 S starts from a high 
value until it suddenly drops to zero at the critical density 
naf Pa 0.00078, cf. bottom part of Fig.g This behavior per- 
sists up to zero temperature, cf. Fig. [6] Vanishing of quantum 
coherence upon crystallization is a general feature in this sys- 
tem and indicates that there is no supersolid phase of indirect 
excitons. If the temperature is above the left triple point the 
superfluid fraction is exactly zero, and compression leads to a 
phase transition from the normal fluid to the crystal phase, cf. 
7s for T = 0.001 Ha* and the change of the PDF in Fig.|4ji-c. 
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Figure 7: Boundaries of the exciton crystal for different dipole mo- 
ments d. Lower abscissa: density range given by D c = 17 and 
r\ = 9.4. Upper abscissa: maximum temperature estimated from 
Id?*. No solid phase exists for d < d c w 9.1 ti£. 



D. Reentrant quantum melting 

Interestingly, if the density is increased further, the exciton 
crystal melts, cf. Fig.|4^,f, this time accompanied by a jump of 
the superfluid fraction from zero to about 0.9. This indicates 
isothermal quantum melting to a (partially) superfluid exciton 
liquid. This occurs at a density of n c a B 2 = 0.0036 ± 0.0003 
corresponding to = 9.4 ± 0.3 and, again, persist to zero 
temperature. At temperatures above the right triple point, 
k^T > 0.001 Ha*, melting and onset of superfluidity are de- 
coupled: first the crystal melts into a normal fluid which be- 
comes a superfluid only at a higher density, cf. Fig. [6] 

Thus the most striking feature of the exciton phase dia- 
gram is the existence of two quantum freezing (melting) tran- 
sitions, even in the ground state. At low-density excitons un- 
dergo pressure crystallization which is characteristic for the 
behavior of dipole systems or, more generally, for neutral mat- 
ter composed of atoms or molecules. In addition, at higher 
densities, there is a second transition: quantum melting by 
compression. While such an effect is absent in conventional 
neutral matter it is ubiquitous in Coulomb systems, including 
the Wigner crystal of the strongly correlated electron gas, ion 
crystals in the core of white dwarf stars and nuclear matter in 
the crust of neutron stars. The existence of this quantum melt- 
ing transition in indirect excitons is due to the peculiar shape 
of the effective potential Vxx- one readily confirms in Fig. [T] 
that at the critical density where the mean exciton-exciton dis- 
tance equals 9.4 a B , Vxx essentially follows the Coulomb re- 
pulsion of the holes (red dashed curve). 



V. CONCLUSIONS 

We have shown that a bosonic many-particle system pos- 
sesses, besides its weakly nonideal Bose condensed gas and its 
superfluid liquid phases also a strongly correlated solid phase. 
Indirect excitons in semiconductor quantum wells have been 
found a favorable candidate due to their long-range pair in- 
teraction and the possibility to achieve strong nonideality by 



controlling the dipole moment with an external electric field. 
Based on first principle PIMC simulations we have computed 
the complete phase diagram in the region of the exciton crys- 
tal. (Quasi-)Long range crystalline order and macroscopic 
quantum coherence are found to be incompatible in an exciton 
crystal - there is no supersolid phase, as long as the crystal is 
free of defects. 



A. Experimental realization 

The results presented above were computed for d = 
13.3 a B . Using values from Refp^ this dipole moment can 
be achieved in a ZnSe quantum well of width L ss 50 nm or 
a GaAs quantum well with L s» 148 nm, both at an electric 
field strength of E = 20 kV/cm. The density interval for the 
exciton crystal is estimated as 1.3- 10 9 cm~ 2 . . . 3.6- 10 9 cm -2 
for GaAs and 8.2 • 10 9 cm~ 2 ... 3.8 • 10 10 cm~ 2 for ZnSe. An 
estimate for the maximum temperature where the crystal can 
exist is obtained from the classical dipole melting curve, 

fc B T dip = c^ (««B 2 ) 3/2 Ha* , (26) 

where c « 0.09,^1 and the critical density n c a B 2 = 0.0036 
is being used. Taking into account that this value is approx- 
imately a factor 2 too high, cf. Fig. [6] we obtain the es- 
timates fc B T max = 0.17K (GaAs) and fc B T max = 0.78 K 
(ZnSe). These parameters are well within reach of current 
experiments. A particular advantage is that the upper density 
limit for exciton crystallization is a factor 16 higher than the 
threshold for an electron Wigner crystal (r s sa 37). A suitable 
diagnostics for the excitonic crystalline phase can be Bragg 
scatteringP21 

B. Dependence of the quantum well width 

Let us now analyze the dependence of the phase diagram 
on the dipole moment d. In semiconductor quantum wells 
the dipole moment can be varied in a broad range by varying 
the QW width or/and the electric field strength. As shown in 
Fig. [7] an increase of d reduces the lower density limit of the 
crystal phase whereas the upper boundary remains unchanged. 
Thus, the crystal phase expands with d, the maximum temper- 
ature T max grows quadratically, cf. Eq. ( |26| > and Fig. [7] Fi- 
nally, there exists a minimum value d c = 9.1 a B where the 
two limiting densities converge, and the exciton crystal phase 
vanishes. 



C. Outlook 

Let us now briefly discuss effects which have been ne- 
glected by the present model, most importantly, disorder and 
thermal relaxation. 

To reduce the effect of the exciton localization at surface 
imperfections we considered the model of a single wide QW 
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(L > 400 A). This allows us to completely neglect the effect 
of 1 monolayer well width fluctuations on the exciton binding 
energy and localization. Some quantitative analysis can be 
found in A. Filinov et al^ In our case, the in-plane size of 
the exciton wavefunction is comparable to the dipole moment 
d = 13.3ctg w 400 Aand is, therefore, of the order of the 
characteristic lateral size of the interface fluctuations ~ 400 A 
(see D. Gammon et alF$. Hence, once the exciton is on the 
top of the defect, the corresponding potential gets significantly 
smoothed. 

In many optical experiments excitons are created in a highly 
non-equilibrium state with a possible coherence and coupling 
to the laser field. Such conditions, certainly, complicate both 
the interpretation of the experiment and the theoretical de- 
scription, and have been studied in detail for polaritons. In 
contrast, we consider an experimental realization, where the 
excitons are created by an optical pulse, which is switched off 
after a short duration, or is periodically repeated with a delay 
of several microseconds, sufficient for the exciton equilibra- 
tion. Fast exciton recombination is prevented by the spatial 
e-h separation due a constantly applied electric field. This sit- 
uation is experimentally feasible as was shown by Z. Vortis et 

aim 



Finally, the most striking feature of the crystal of indirect 
excitons, confirmed by the simulations, is two quantum melt- 
ing transitions which persist at zero temperature: at low den- 
sities it melts by expansion whereas at high densities it melts 
when being compressed. The origin of this unusual and rich 
phase diagram has been traced to the non-trivial form of the 
exciton interaction potential. With it the exciton solid com- 
bines features of conventional neutral matter (exhibiting crys- 
tallization by compression) and Coulomb matter (quantum 
melting by compression), as found for instance in exotic com- 
pact stars. 
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